i_tmp=int32((x(i)-minx)/xfd)+1; j_tmp=int32((y(i)-miny)/yfd)+1; z_tmp=int32((z(i)-minz)/zfd)+1; cntxy=tjxy(i_tmp,j_tmp ); cntxz=tjxz(i_tmp,z_tmp ); cntyz=tjyz(j_tmp,z_tmp ); cntxyz=tjxyz(i_tmp,j_tmp,z_tmp ); tjxy